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polarization survey: implications for CMB observations 
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ABSTRACT 



Aims. The aim of the present analysis is to improve the knowledge of the statistical properties of the Galactic diffuse synchrotron emission, 
which constrains sensitive CMB anisotropy measurements. 
£^ Methods. We have analysed the new DRAO 1.4 GHz polarization survey together with the Stockert 1.4 GHz total intensity survey and derived 
\^ the angular power spectra (APSs) of the total intensity, the polarized emission, and their cross-correlation for the entire surveys and for three 
£^ , low-intensity regions. 

Results. The APSs of the diffuse synchrotron emission are modelled by power laws. For the E and B modes, a slope of a ~ [-3.0, -2.5] for the 
multipole range ~ [30, 300] is found. By the extrapolation of these results to 70 GHz, we can estimate the Galactic synchrotron contamination 
of CMB anisotropics, and we find results that are compatible with the ones coming from WMAP 3-yr data. In the low-intensity regions, 
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O 1 the cosmological primordial B mode peak at I ~ 100 should be clearly observable for a tensor- to- scalar ratio T/S > 0.5 and a synchrotron 
^ temperature spectral index J3 ~ -3. Its detection is also possible for T/S ~ 0.005 and J3 ~ -3, in case a separation of the foreground from the 
^ j CMB signal could be achieved with an accuracy of ~ 5 - 10%. For the TE mode, a mask excluding \b ga i\ < 5° (for p ~ -3.0) or \b ga i\ < 20° 
• • . (for p ~ -2.8) from the surveys is sufficient to render the foreground contamination negligible, thus confirming the ability of WMAP to have a 
^ clear view of the temperature-polarization correlation peak and antipeak series. 

CZ Key words. Polarization - Galaxy: general - Cosmology: cosmic microwave background - Methods: data analysis. 



1. Introduction 

The polarized Galactic diffuse synchrotron radiation is ex- 
pected to be the major foreground at v < 70 GHz on inter- 
mediate and large angular scales (0 > 30') at medium and 
high Galactic latitudes. At about 1 GHz, the synchrotron 
emission is the most important radiative mechanism out of 
the Galactic plane. Consequently, radio frequencies are the 
natural range for studying it, though it might be affected 
by Faraday rotation and depolarization. Before the DRAO 
survey became available, the Leiden polarization surveys 
( |Brouw & Spoelstr a 1976| ) provided a sky coverage, allowing 
to estimate the Galactic synchrotron APS on intermediate and 
large scales ( |La Porta & Bur igana 2006 }. However, these sur- 
veys have sparse sampling, low sensitivity, and a good signal- 
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to-noise ratio only for the brightest regions in the sky. A 
new linear polarization survey of the northern celestial hemi- 
sphere at 1.42 GHz with an angular resolution Ohpbw - 36' 
has recently been completed using the DRAO 26 m tele- 
scope (Wolleb en et al. 20 06). The survey has a spacing, S , of 
15' in right ascension and from 0?25 to 2°. 5 in declination, 
which requires interpolation to construct equidistant cylindri- 
cal (ECP) maps with 6 pixd = 15' (Wolleben 20051) . The fi- 
nal map has an rms-noise of about 12 mK, which is unique 
so far in terms of sensitivity. The polarized intensity map ap- 
pears cold and patchy at high Galactic latitudes, except for 
the "North Polar Spur" (NPS), extending far out of the plane 
at about l ga \ ~ 30° 1 . We used the polarization data in com- 
bination with the Stockert total intensity map at 1.42 GHz 
(Reich 1982; Reich & Reich 1986), having the same angu- 
lar resolution, similar sensitivity, and Q s ^ 15', to investigate 



1 The NPS has been extensively studied (see Egger & 
Aschenbach (1995) and references therein) and is interpreted as 
the shock front of an evolved local supernova remnant. 
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Fig. 1. Gnomonic projection of the cold patches representative 
of the diffuse Galactic synchrotron emission at 1.4 GHz (total 
intensity - top panels; polarization - bottom panels). 



the statistical properties of the Galactic diffuse synchrotron 
emission. The major purpose of the analysis is the compar- 
ison of the foreground fluctuation properties with those of 
the CMB anisotropy field, as drawn in standard cosmologies. 
Therefore, we described the anisotropy field statistical prop- 
erties in terms of its angular power spectrum (APS) 2 (see 
IPeebles 19931 IKamionkowski et al. 19971 l^ldarriag alUoT] ) 
and adopted HEALPix 3 (Gorski et al. 2005| as the tessellation 
scheme for the sphere. A detailed description of the algorithm 
to convert the original ECP maps into HEALPix maps was given 
by La Porta et al. (2005). 

2. Selected areas and data analysis 

We selected three areas with low polarized intensities (see 
Fig. [I}, which are identified by the following Galactic coor- 
dinates: Patch A - 180° < £ gaX < 276°, b gal > 45°; Patch B - 
193° < £ gal < 228°, b ga i < -45°; Patch C - 65° < £ gal < 180°, 
b ga i ^ 45°. Being interested in the diffuse Galactic synchrotron 
emission at 1.4 GHz, we subtracted discrete sources (DSs) from 
the total intensity map (very few sources are visible in polariza- 
tion out of the plane). We performed a 2-dimensional Gaussian 
fitting using the NOD2-software library (Hasla m~1974l» . which 
also allowed us to estimate the diffuse background. We are 
confident that at \b ga i\ ~ 45° DSs with peak flux densities ex- 
ceeding ~ 1 Jy have been removed that way (see Burigana et 
al. 2006 for a map of the subtracted DSs). 
We then derived the APSs of the fluctuation fields (i.e., after 
subtraction of the mean value from the maps) for the three se- 
lected areas and the whole survey coverage using the HEALPix 
facility anafast. Figure |2| shows the APSs obtained for the 
selected areas in both temperature and polarization. 

In all cases, the polarization APSs are rather similar 4 . 
This fact might indicate that depolarization due to differential 

2 The angular scale 6 and the angular power spectrum multipole £ 
are related by £ ~ 18O°/0[°]. 

3 http://healpix.jpl.nasa.gov/ 

4 In particular, Cf - Cf We then introduce the quantity Cf ,B = 
0.5 -(Cf + Cf , which will be used in the following best-fit analysis. On 
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Fig. 2. APSs of the selected cold regions (solid lines — > Cj 7 ; 
dashes —> Cf ; dots —> Cf\ three dots-dashes — > Cf). 
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Fig. 3. APSs of the selected cold regions and correspond- 
ing best-fit curves (dot-dashed lines). From the bottom: C 77 , 
20 • Cf and 50 • Cf B . In the case of Cj 7 , cf^V^ and 
c src e -{cr b £f are a i SQ displayed. 
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Table 1. Least-square best-fit parameters describing the APSs 
of the selected patches in the range 30 < £ < 300. The errors 
are given so that the upper (resp. lower) values correspond to 
the best-fit parameters set with the flattest (resp. steepest) slope. 



Faraday rotation should not be relevant in such sky regions at 
the investigated angular scales. In fact, rotation measure (RM) 
maps ( Joh nston-Hollitt et al. 20041 IDineen & Coles 20051 ob- 
tained interpolating RM data of extragalactic sources show 
very low values in correspondence to such areas. However, 
the degree of polarization is on average a few percent, 
well below the maximum theoretical value of ~ 75% 
(Ginzburg & Syrovatskii 1965| ). The reason for the low frac- 
tional polarization is not clear. One possibility is that de- 

the contrary, the E and B modes of the primordial CMB anisotropics 
differ largely as they are induced by different mechanisms (e.g., Seljak 
& Zaldarriaga 1997; Kosowsky 1999). 
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polarization other than differential depolarization is present 
(e.g.. lSokoloff etal. 19981) . 

The maps show the Galactic synchrotron emission and non- 
subtracted DSs convolved with the telescope beam and con- 
taminated by noise. We then fitted the APSs of both total inten- 
sity and polarization (namely Cj 7 , C^ 1 , and C^ ,B ) by modelling 
them as Q = (Cf nch + c src ) • + c WN . We exploited the power 
law approximation cf nch = k • £ a and assumed a symmetric 
Gaussian beam, i.e., a window function Wf = e~ (cr ^ 2 , where 
o~b - 6hpbw(J&&)I V81n2 . The contribution of non- subtracted 
DSs has been simply modelled with a constant term accord- 
ing to the formalism of Poisson fluctuations from extragalactic 
point sources (Toffolatti et al. 19981 . expected to represent the 
largest fraction of the non- subtracted DSs. From extragalactic 
source counts at 1.4 GHz (Pran doni et al. 200Tb we estimate 
c src ~ o.09 mK 2 for Cj 7 (adding or subtracting quoted Icr er- 
rors to the best-fit number counts we (conservatively) get c src ^ 
[0.05 - 0.3] mK 2 ). A guess of the noise contribution is given 
by Cj N ~ An • cr p i X /N P i X , where N P i X is the number of pixels 
in the HEALPix map and cr pix = cr pixMkL?ix ~ N~ l/2 • o> x , E cp, 
N being the number of the ECP pixels corresponding to each 
HEALPix pixel at a given latitude (for the total intensity maps 
of the selected areas, ct^healpix ~ 13 mK corresponding to 
C WN ~ 0.003 mK 2 ). The best-fit results for all patches are 
listed in Table [T] and shown in Fig. [3] For temperature, the 
parameter range for the Galactic synchrotron emission APS is 
a ~ [-2.71,-2.60] for the slope and log 10 * ~ [4.067,4.339] 
(k in mK 2 ) for the amplitude. The DS contribution is in the 
range c src ~ [0.023,0.047] mK 2 , which is consistent with the 
above source counts estimate, within the errors 5 . For polariza- 
tion, the derived slopes are generally steeper, varying in the 
interval a ~ [-3.02, -2.62] for Cf and a ~ [-3.05, -2.55] for 
Cf ,B . The DS contribution is on average much lower than in 
temperature, since it is compatible with zero 6 . 



3. Discussion: implications for CMB observations 

We have extrapolated the recovered APS of the E and B modes 
to 70 GHz and compared it with the APS of the CMB polariza- 



5 The bulk of the factor ~ 2 of discrepancy between the value of c src 
recovered by our fit and that predicted from best-fit number counts can 
in fact be produced by the survey sky sampling (0 S ) of 15'. For exam- 
ple, considering a lower limit corresponding to ^ 15'/ V2 (instead of 
0) in the integral over if/ in Eq. (2) of Toffolatti et al. (1998), c src de- 
creases by a factor 1.53. 

6 The best-fit results may suggest a polarization degree (obtained 
considering the contribution of the subtracted DSs, ~ 0.05 - 0.2 mK 2 , 
to the temperature APS also) considerably higher than ~ 2%, the 
value found for NVSS extragalactic (mainly steep spectrum) sources 
(Mesa et al. 2002; Tucci et al. 2004). It may imply a presence (or 
a combination) of spurious instrumental polarization at small scales, 
of a significant contribution from highly polarized Galactic sources 
^Manchester et al. 19981 non- subtracted in the maps, or of a flattening 
of the diffuse synchrotron polarized emission APS at i > 200 - 250 in 
higher resolution data on smaller sky areas (Baccigalupi et al. 2001; 
Carretti et al. 2006). 



tion anisotropy 7 for a ACDM model including scalar and tensor 
perturbations compatible with the recent WMAP 3-yr results 8 
(Spergel et al. 2006). The frequency range between 60 GHz and 
80 GHz seems the less contaminated by Galactic synchrotron 
and dust foregrounds in both temperature (Bennett et al. 2003) 
and polarization (Page et al. 2006) at angular scales > 1°. The 
accurate measure of the E mode is of particular relevance for 
breaking the existing degeneracy in cosmological parameter es- 
timation, when only temperature anisotropy data are available 
(e.g., Bond et al. 1995, Efstathiou & Bond 1999). The detection 
of the primordial B mode is of fundamental importance for test- 
ing the existence of a stochastic cosmological background of 
gravitational waves (e.g., Knox & Turner 1994). The results of 
our comparison have been displayed in Fig.|4|for two choices of 
the temperature spectral index (J3 = -2.8,-3). The APS extrap- 
olated from the entire DRAO survey is also shown. It exceeds 
the CMB E mode even at the peak at £ ~ 100 for p = -2.8. For 
J3 = -3, the two APSs are almost comparable. Figure @] also 
shows the APS extrapolated from the DRAO survey excluding 
the region \b ga i\ < 20° and adopting J3 = -3 9 . Such a sky mask 
reduces the Galactic APS below the CMB E mode for £ > 50. 
In this case, the CMB B mode peak at £ ~ 100 is comparable to 
(or exceeds) the synchrotron signal for tensor- to- scalar ratios 
T IS > 0.5. For the three patches and T/S > 0.5, the cosmo- 
logical B mode peak at £ ~ 100 is comparable to (or exceeds) 
the Galactic synchrotron signal extrapolated with j3 -2.8, 
while it is larger by a factor > 2 (in terms of VQ) f° r J3 - -3. 
Furthermore, a separation of the Galactic synchrotron polar- 
ized signal from the CMB one with a ~ 5 - 10% accuracy (in 
terms of VQ) would allow to detect the CMB B mode peak 
at £ ~ 100 even for T/S ~ 0.005 if p -3. Similar results 
for the detection of the B mode peak at £ ~ 100 have been 
inferred from 1 .4 GHz polarization observations of a small re- 
gion with exceptionally low Galactic foreground contamina- 
tion ICarretti et al. 2 006). though at £ ~ few x 100 the CMB B 
mode is expected to be dominated by the B mode generated by 
lensing (Zaldarriaga & Seljak 1998). Finally, we note that for a 
sky coverage comparable with those of the considered patches 
(~ 3%), the cosmic and sampling variance does not signifi- 
cantly limit the accuracy of the CMB mode recovery at £ larger 
than some tens. 

The CMB TE correlation constrains the reionization his- 
tory from the power at low multipoles and the nature of primor- 
dial fluctuations from the series of peaks and antipeaks (e.g., 
Kogut 2003; Page et al. 2006). In Fig. we compare the TE 
mode APS of the Galactic emission extrapolated to 70 GHz 
with the CMB one. The APS of the whole DRAO survey indi- 
cates a significant Galactic contamination at £ < fewx 10, even 
for/? -3. The use of suitable Galactic masks (e.g., excluding 
regions at \b gat \ < 5° for /3 ~ -3 or \b gat \ < 20° for /3 ~ -2.8) 

7 We used the CMBFAST code (version 4.5.1) for the computation of 
the CMB APS (http://www.cmbfast.org/). 

8 http://lambda.gsfc.nasa.gov/ 

9 This APS is consistent with the WMAP foreground polarization 
APS at 63 GHz (see Fig. 17 in Page et al. 2006), thus supporting an 
effective 1.4-70 GHz polarization spectral index J3 -3 far out of 
the Galactic plane. 
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Fig. 4. Comparison between the E and B modes of the CMB 
anisotropy at 70 GHz and of the Galactic polarized synchrotron 
foreground. Dotted and solid lines: CMB E and B mode for two 
different tensor- to-scalar ratios (thin line T/S = 0.005, thick 
line T/S = 0.5). Dashed lines: uncertainty associated with the 
cosmic and sampling variance for a sky coverage of ~ 3% and 
a binning of 10% in £. Three dots-dashes: average of the E 
and B modes extrapolated from the whole DRAO survey for 
spectral indices f3 - -2.8 and -3, respectively. The anafast 
results have been divided by the window function to correct for 
beam smoothing. Thick long dashed line: the same as above, 
but with a mask at \b ga i\ < 20 for/3 = -3. Dashed (dot-dashed) 
power law: best-fit result corresponding to Patch B (resp. C), 
rescaled in frequency as above. The results are here in terms of 
thermodynamic temperature. 

greatly reduces the Galactic foreground contamination even at 
low multipoles. For the three cold patches, the Galactic TE 
mode is negligible compared to the CMB one, independently 
of the adopted J3. The TE mode antipeak at I ~ 150 turns out 
to be very weakly affected by Galactic synchrotron contamina- 
tion in all cases. 
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